CERN-TH/98-151 



Q2-EV0LUTI0N OF SPIN-DEPENDENT 
PARTON DENSITIES 
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Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

We discuss the NLO evolution of quark transversity densities and of 
the parton distribution function for linearly polarized gluons in a linearly 
polarized hadron. A supersymmetric relation between the NLO evolution 
kernels for transversity and for linear polarization is found. We also study 
the implications of NLO evolution for Soffer's inequality and the prospects 
of measuring transversity densities in polarized Drell-Yan at RHIC. 



1. Introduction 

Experimentally, the vast majority of data sensitive to parton densities 
have been taken without fixing the polarization of the initial beams or the 
target. The densities extracted in this way are usually refered to as the 
'unpolarized' parton distributions f{x,Q'^) (/ = q,q,g)- 

Within roughly the last decade, also more and more data have be- 
come available that are sensitive to the 'longitudinally' polarized ('helicity- 
weighted') parton densities of the nucleon. The tool to obtain such informa- 
tion has (almost exclusively) been deep-inelastic scattering (DIS) of longi- 
tudinally polarized leptons and nucleons. The spin-asymmetry measured in 
such reactions is related to the probability for finding a certain parton-type 
with positive helicity in a nucleon of positive helicity minus the probability 
for finding it with negative helicity. These densities, denoted as Af{x,Q'^) 
if = QyQyd)^ contain information different from that contained in the more 
familiar unpolarized ones. 

For a transversely polarized spin-1/2 hadron one can define a further 
quark density |, i in very much the same way as the longitudinally po- 
larized quark distributions, by taking differences of probabilities for finding 
quarks with transverse spin aligned and anti-aligned with the transverse 
hadron spin. These densities are called 'transversity' densities and are de- 
noted by At/(x, Q^). It turns out that - unlike the situation for unpolarized 
and longitudinally polarized densities - there is no gluonic analogue of quark 
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transversity [0. This is due to angular momentum conservation: transver- 
sity densities are related to helicity-flip amplitudes. A gluonic helicity-flip 
amplitude would require the hadron to absorb two units of helicity, which 
no spin-1/2 target can do. 

The transversity distributions are completely unmeasured so far since 
they cannot be directly accessed in DIS. It seems certain, however, that 
measurements of the Aj^/ will be attempted at the future polarized proton- 
proton collider RHIC at Brookhaven The most suitable candidate for 
such measurements is believed to be Drell-Yan dimuon production 

The Ay/ complete the twist-2 sector of parton densities of spin-1/2 
hadrons. Nevertheless, we have not yet depleted the full set of parton den- 
sities that can be defined if spin is taken into account. There finally is a 
further spin-dependent gluon distribution that to some extent can be re- 
garded as the gluonic counterpart of quark transversity. Unlike the 'helicity' 
density Ag that describes circular polarization of the gluon, it is encoun- 
tered if the gluon is linearly polarized [§, |5|, ^. The density is denoted by 
ALg{x,Q'^) and exists only in a linearly polarized hadron (or photon Q), 
which therefore has to have spin > 1. There is no quark distribution in this 
case m. Even though a measurement of A^g does not seem very realistic at 
the moment, it possesses some interesting theoretical aspects which justify 
its analysis. 

Table 1 summarizes the parton densities we have defined. 

It is important to realize that each set of parton densities in Tab. 1 (i.e.. 



polarization 


quarks 


gluons 


unpolarized 


q = ql + q+=q]+q\ 




long, polarized 


Aq = qX- q+ 


Afi- = si - 9+ 


transversity 


^rq = q] -q\ 




linearly pol. glue 




^L9 = 9l-4 



Table 1. List of twist-2 quark and gluon densities including spin-dependence. We have 
suppressed the ubiquitous argument (a;, Q^) of the densities. Note that 'q always runs 
over quarks as well as over antiquarks. Labels +, — denote helicities, f, j transverse 
polarizations, and x (y) stands for linear polarization along the x (y) axis, where the 
particle is moving along the z-direction. Subscripts refer to partons and superscripts 
to the parent hadron. 
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each of the rows of Tab. 1) is subject to its own set of evolution equations. 
For instance, the evolution of the longitudinally polarized densities proceeds 
independently from that of the unpolarized partons, and so forth. In this 
way, we are led to introducing separate sets of evolution kernels (splitting 
functions) for each type of polarization. 

The Q^-evolution of the unpolarized densities has been worked out up 
to NLO accuracy of QCD already a long time ago 0, |^, and it has be- 
come standard since about ten years to analyse the unpolarized data within 
the NLO framework. The LO evolution of the Af{x,Q^), the ATf{x,Q'^), 
and of ALg{x,Q^) has also been known for a long time 0, |lO|, ^, whereas 
the derivation of the NLO evolution kernels for longitudinally polarized par- 



tons has been a more recent development [11, 12[. Very recently, the NLO 
splitting functions for transversity were derived within three independent 
calculations ||l^, |l^]. In this paper, we will discuss the NLO evolution of 
the At/, and its implications for an inequality between the /, Af and 
Axf derived by Soffer |jl^. We will also for the first time present the NLO 
evolution kernel for Aig- 



2. NLO evolution equations 

Since there are no gluons involved in the case of the transversity distribu- 
tions, their evolution equations reduce to simple non-singlet type equations. 
Introducing 

ArqliQ') ^ Arq^iQ^) ± At^ (Q') , (1) 
where the Mellin moments are defined by 



Jo 



(2) 



one has the evolution equations (see, for example, [ p^ ) 
d 



di^Q2^Tqim = ATP^^^^{as{Q'))ATql{Q') (3) 

for all flavours. The splitting functions AxPgg -i-{as{Q'^)) are taken to have 
the following perturbative expansion: 

ATP,"„i(a.) = (I) A^P(°)- + (g)' A^PW- + . . . . (4) 
As indicated, ATPqq,+ and AxPqq- are equal at LO. 
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Similarly to (^) one has for linearly polarized gluons: 
d 



dlnQ2 



Ai5"(Q2) = Az.P!l(a,(Q2))Ai5"(Q2) , (5) 



where A^P^ has a peturbative expansion analogous to @. These NLO 



evolution equations can be easily solved; for details see, for instance, [13|. 



3. Calculation of splitting functions 

Historically, the evolution kernels for parton densities have been cal- 
culated within two rather different methods. On the one hand, the very 
powerful, yet quite formal, technique of the 'Operator-Product-Expansion' 
(OPE) has been applied. Here the evolution kernels are derived as anoma- 
lous dimensions of matrix elements of local operators. The other method 
is more intuitive and relies on parton model ideas and on the factorization 



of mass singularities in a physical (ghost-free) gauge |17, 0. This is the 
method we will use. Very roughly, the strategy goes as follows: The NLO 
splitting functions are related to the residues of the collinear singularities in 
the next-to-NLO (NNLO) calculation of a partonic subprocess cross section 
for a 'physical' process. Being the pole part of this cross section, the singular 
terms are manifestly gauge-independent. This means that we can employ 
a gauge in their calculation that simplifies the extraction of the mass sin- 



gularities as much as possible. As was shown in [17|, working in a physical 
gauge, like the light-cone gauge, reduces the number of singular graphs at 
a given order significantly. Thus, a full NNLO calculation is clearly not re- 
quired. In particular, only ladder-like diagrams, corresponding to a parton 
cascade, contribute to the pole part in such gauges, so that the cross section 
a for any physical partonic process can be diagrammatically expanded into 
a sum of ladders of 'two-particle-irreducible' kernels that are individually 
finite. Mass singularities only occur on the lines connecting the kernels (i.e., 
on the 'sides' of the ladders) and can therefore be projected out easily using 
a projector V that converts these lines into on-shell physical states. Thus, 
the factorization of mass singularities can be achieved: the cross section a 
decomposes into a finite part ap, and a function T that contains all the 
collinear singularities and is universal since no process-dependence is left 
over in it. The function T is then directly related to the splitting functions 
one is looking for [0, ^ . 

For the case of transversity, the projector AtV is found to be 

^tV ~ , (6) 
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where k is the momentum of the particle emerging from the top of a kernel, 
n is the gauge vector with n • A = 0, = for the light-cone gauge, 
and s is the transverse spin vector. More details on the calculation in the 
transversity case can be found in [p!^ ]. 
For linearly polarized gluons one has 

where the polarization vectors 

= (0,1,0,..., 0) , ej; = (0,0,1,0,...,0) (8) 

have non-vanishing entries only in the x and y components, respectively. A 
more detailed account of our calculation of A^P^^^ will be given in a future 



publication [18|. 



4. Results 



Our results for ATPql^±{x) can be found in |13] and need not be repeated 



qq 



here. The one for A^Pgo is new and reads in the MS scheme 



N, 



c 



fiT 1 — 

— + - ln2 X - 2 In X ln(l - x) - y j 6lP^^'^ (x) 
+ S2ix)5LP^fi-x) + (^^ + 3C(3)^ 6{1 



18 ' 2 
1 — x^ 



where 



+ NcTf 



6x 

2— + .(1-.) 



2x 



(9) 



(10) 



(l-x)+ ' 

1+^ , ,1 — z . 
— ln( ) 

^ z z 

l + x 

-2Li2(-x) - 21nxln(l +x) + - In^ x , (11) 

2 6 
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with Li2(x) being the dilogarithm. Furthermore, we have as usual Cp = 4/3, 
Nc = 3, and Tf = nf/2, nj being the number of active flavours. Finally, 
C(3) 1.202057. We note that in |] a result for the part ~ CpTf in (|) at 
X < 1 was presented, which corresponds to the two-loop splitting function 
for transitions between linearly polarized gluons and photons. The result of 

disagrees with ours by a an additional factor 1/x in |^. The calculation 
of Q therefore has to be in error, since a behaviour ~ 1/x^ of the splitting 
function - more singular than the unpolarized one - cannot be correct. 

There are two aspects of the result (^) that deserve further attention. 
Firstly, the small-x behaviour of the splitting function for linear polariza- 
tion changes quite dramatically when going from LO to NLO. At LO, the 
splitting function behaves ~ x as x ^ 0. At NLO, there are terms ~ 1/x in 
the splitting function (as one also encounters in the unpolarized case); we 
have 



AiP«(x)p.;l 
"'^ 6x 



+ 0{x) 



[X 



0) . (12) 



We note that all logarithmic terms ~ xln^ x cancel out in this limit. 

The other interesting point concerns the 'supersymmetric' limit Cp = 
Nc = 2Tf = N |l^ which has already been investigated for the unpolar- 
ized and the longitudinally polarized NLO splitting functions in 11, 12| 



and for the 'time-like' ones in |21|. One first of all finds that in the super- 
symmetric limit the LO splitting functions for quark transversity and for 
linearly polarized gluons become equal [P, ^: 



ATPr( 



X] 



N 



2x 3,,^ 



(13) 



Thus, we have found the supersymmetric 'counterpart' of transversity: lin- 
ear polarization (see also [|l9|). This nicely completes the supersymmetry 
relations found in the case of unpolarized and longitudinally polarized par- 
ton densities and fragmentation functions. The fact that two rather differ- 
ent polarizations are linked here is not too surprising as both of them are 
transverse in a certain sense. 

To see whether the supersymmetric relation also holds at NLO, we have 
to transform the splitting functions to a regularization scheme that respects 
supersymmetry, dimensional reduction. The procedure here follows closely 
that discussed in |12, 21 1. The task is simplified by the fact that there 
are no parts ~ e in the d-dimensional LO splitting functions for transver- 
sity or linearly polarized gluons at x < 1. Such terms - which are always 
absent in dimensional reduction, but can be present for dimensional regu- 
larization ~ are the reason why, for instance, the longitudinally polarized 
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NLO splitting functions do not directly satisfy the supersymmetric relation 
in dimensional regularization, but have to be transformed to dimensional 



reduction first (see for a more detailed discussion). The fact that the 



LO e-terms are absent at x < 1 in our case, means that our final results for 
the NLO transversity splitting function AxPqg^ o f J3[ and for the NLO 
splitting function for linearly polarized gluons, Eq. (|^7already coincide for 
X < 1 with their respective MS expressions in dimensional reduction. We 
can therefore immediately compare these expressions in the supersymmetric 
limit. Indeed, we find for Cp = Nc = 2Tf = N: 

ATP^llix) ^ AlP^I\x) (x<1). (14) 

The satisfaction of the supersymmetric relation at x = 1 is trivial; see [|l2| , 
where the appropriate factorization scheme transformation to dimensional 
reduction at x = 1 is given. 

We thus find that indeed the supersymmetric relation between the trans- 
versity splitting function and the evolution kernel for linearly polarized glu- 
ons is also verified at NLO. We note that the relation involves as expected 
the ' -I— combination' ATPqg\ = A^Pgg^ + Aj-Pg^^ which corresponds to 'sin- 
glet' evolution. Finally, also note the miraculous cancellation of the terms 
~ 1/x in Aj^Pgg^ in the supersymmetric limit (see Eq. (|12|)) - after all, such 
terms are not present in AxPgg^^^- 



5. SofFer's inequality and the transversely polarized Drell Yan 

process 

The unpolarized, longitudinally and transversely polarized quark dis- 
tributions {q, Aq, Axq) of the nucleon are expected to obey the rather 
interesting relation 

2\ATqix)\ < q{x) + Aq{x) (15) 



derived by Soffer |15|. It has recently been clarified that Soffer's inequation 
is preserved by leading order (LO) QCD evolution, i.e. if (|l5|) is valid at 
some scale Qq, it will also be valid at Q > Qo To NLO the situation is 
not as simple. The parton distributions are now subject to the choice of the 
factorization scheme which one may fix independently for q, Aq and A^g. 
One can therefore always find 'sufficiently incompatible' schemes in which 
a violation of ([l^) occurs. However, it was shown in fl^ with analytical 
methods that the inequality for valence densities is preserved by NLO QCD 
evolution in a certain 'Drell- Yan scheme' in which the NLO cross sections for 
dimuon production maintain their LO forms, and also in the MS scheme. In 
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this section, which is taken from |23|, we shall show numerically that Soffer's 
inequation for sea quarks is also preserved under NLO (MS) evolution. We 
note that there is also a recent analytical proof of the preservation of Soffer's 



inequality IJ at NLO. 

For our study, we will assume saturation of Soffer's inequality at the 
input scale for parton evolution. Our choice for the r.h.s. of ([l5|) w ill then 
be the NLO MS radiative parton model inputs for q{x, Qq) of |25|| and for 
the longitudinally polarized Aq^x^Q^) of the 'standard' scenario of |26] at 



Qo = tJ^'jsiLO ~ ^-^^ GeV^. For simplicity we will slightly deviate from the 
actual Qq) of in so far as we will neglect the breaking of SU(2) 
in the input sea quark distributions originally present in this set. This 
seems reasonable as SU(2) symmetry was also assumed for the Ag(x, Qq) 
of [|2^, which in that case was due to the fact that in the longitudinally 
polarized case there are no data yet that could discriminate between An 
and Ad. We therefore prefer to assume Atu{x,Q'^) = Ar(i(x,Qo) ^-Iso for 
the transversity input. 

In order to numerically check the preservation of (p]5|). Fig. [l| shows the 
ratio 

as a function of x for several different values for q = = u^, u = 
{u+ - U-)/2, d., = d-, d = {d+ - d-)/2 (cf. Eq. (D). If NLO evolution 
preserves Soffer's inequality, then Rq{x,Q'^) should not become larger than 
1 for any > Qq. As we already know from [^] this is the case for 
the two valence distributions. Figure || confirms that Soffer's inequality is 
indeed also preserved for sea distributions. Furthermore, we see in Fig. |l| 
that evolution leads to a strong suppression of Rq{x, Q"^) at small values of 
X, in particular for the sea quarks. 

As is obvious from (|l5|), Soffer's inequality only restricts the absolute 
value of the transversity distribution. Therefore, we are free to choose the 
sign when saturating (|l5|), and we have to check the results for the two 
distinct possibilities /S.Tqv{x,Q'^) > 0, ATg(x,Qo) > ^'^d /S.Tqv{x,Q'^) > 
0, ATq{x, Qq) < 0. Our results do not noticeably depend on the actual 
choice. Neither does it matter whether we decide to saturate Soffer's in- 
equality for the valence densities qv or for the full quark {q = qv + q) distri- 
butions at the input scale. 

One can utilise a saturated Soffer inequality to derive upper bounds 
on the transverse double-spin asymmetry Att to be measured in trans- 
versely polarized Drell-Yan muon pair production at RHIC. Att is defined 
as Att = d5a/da where the polarized and unpolarized hadronic cross sec- 
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Fig. 1. The ratio Rq{x,Q'^) as defined in ( [iq ) for q — Uv,u,dv,d and several fixed 
values oi Q^. 



tions are 



d5a = - 
2 



da 



TT 



da 



n 



d(T^-i da^'' + da^^ 



(17) 



We employ the same input distributions as before, along with the same value 
for the initial scale Qq. We choose Qf = Qr = M for the factorization and 
renormalization scales, where M is the invariant mass of the muon pair. 
Further details of the calculation of the Drell-Yan cross section to LO and 
NLO can be found in |23|. 

Figure ^ shows the transversely polarized pp cross section and the 'max- 
imal' double-spin asymmetry Att at LO and NLO for ^/S = 150 GeV, 
corresponding to the RHIC pp collider. We note that the region 9 GeV 
< M < 11 GeV will presumably not be accessible experimentally since it 
will be dominated by muon pairs from bottomonium decays. The predicted 
maximal asymmetry is of the order of a few per cent. The expected statisti- 
cal error bars have been included in the plot; they have been calculated for 
70% beam polarization, an integrated luminosity £ = 240 pb~^, and 100% 
detection efficiency. One concludes that asymmetries of this size should be 



well measurable at RHIC. Similar asymmetries are found [23| for another 
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4 6 8 10 12 14 16 18 20 6 8 10 12 14 16 18 20 
M [GeV] M [GeV] 

Fig. 2. NLO and LO maximal polarized Drell-Yan cross sections and asymmetries 
for RHIC at = 150 GeV. The error bars have been calculated assuming £ = 
240 pb~^, 70% beam polarization and 100% detection efSciency. 

conceivable experimental situation, namely a proposed fixed-target exper- 
iment at HERA that would utilise the possibly forthcoming polarized 820 
GeV HERA proton beam on a transversely polarized target, resulting in 
V5 = 40 GeV. 

Figure |^ shows similar results for the high-energy end of RHIC, ^/S = 
500 GeV, where the integrated luminosity is expected to be £ = 800 pb~^. 
It turns out that the asymmetries become smaller as compared to the lower 
energies, but thanks to the higher luminosity the error bars become rela- 
tively smaller as well, at least in the region 5 GeV ^ M < 25 GeV where 
the errors are approximately 1/10 of the maximal asymmetry. One can also 
clearly see in Fig. |^ the effect of Z exchange and the Z resonance. 

A comparison of the LO and NLO results in Figs. |2| and ^ answers one 
key question concerning the transversely polarized Drell-Yan process: our 
predictions for the maximal Att show good perturbative stability, i.e. the 
NLO corrections to the cross sections and Att are of moderate size, albeit 
not negligible. 

Finally, we also find a significantly reduced scale dependence of the 
results when going from LO to NLO. This is shown in Fig. ^ for the case 
^fS = 150 GeV. We plot here the maximal asymmetry in LO and NLO, 
varying the renormalization and factorization scales in the range M/2 < 
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Qf = Qr ^ 2M. One can see that already the LO asymmetry is fairly 
stable with respect to scale changes, which is in accordance with the finding 
of generally moderate NLO corrections. The NLO asymmetry even shows 
a significant improvement, so that Att becomes largely insensitive to the 
choice of scale. 

Acknowledgement: I am grateful to O. Martin, A. Schafer and M. Strat- 
mann for a fruitful collaboration. 




Fig. 3. As in Fig. |, but for VS = 500 GeV and C = 800 pb'^ 
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Fig. 4. Scale dependence of the LO and NLO asymmetries at ^/S = 150 GcV. 
The renormaHzation and factorization scales were chosen to be Qn — Qf — 
M/2, M, 2M. The solid Hne is as in Fig. H 
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